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ABSTRACT 

We map out calcium II & sodium I absorption (Fraunhofer H, K & D lines) induced by 
both the interstellar medium and the circumgalactic medium of the Milky Way. Our 
measurements cover more than 9000 deg 2 and make use of about 300, 000 extragalactic 
spectra from the Sloan Digital Sky Survey. We present absorption maps for these two 
species and then compare their distributions to those of neutral hydrogen and dust. 
We show that the abundance of Na I with respect to neutral hydrogen stays roughly 
constant in different environments, while that of Ca II decreases with hydrogen column 
density. Studying how these tracers vary as a function of velocity, we show that, on 
average, the N-n&i/N ca.ii ratio decreases at higher velocity with respect to the local 
standard of rest, similar to the local Routly-Spitzer effect but seen on Galactic scale. 
We show that it is likely caused by higher gas/dust density at lower velocity. Finally, 
we show that Galactic Ca II and Na I absorption needs to be taken into account for 
precision photometry and, more importantly, for photometric redshift estimation with 
star forming galaxies. Our maps of Ca II and Na I absorption are publicly available. 
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1 INTRODUCTION 


Various advances in our understanding of the Universe have 
been enabled by sky surveys mapping emission and absorp¬ 
tion features across the sky. Historically, William Herschel 
pioneered this approach by mapping the density of stars 
and revealing the existence of ’’dark holes” (Herschel, 1785), 
known today as dust clouds. The mapping of globular clus¬ 
ters across the sky allowed Harlow Shapley to infer the disk¬ 
like geometry of the Milky Way about one hundred years 
ago. In the sixties, radio surveys opened up a new win¬ 
dow giving u s access to the distribut ion of hydrogen in the 
Galaxy (e.g., iKalberla fe Kerrill2009l ). Over the past twenty 
years, IR obs ervations have provid ed us with the distribution 
of dust (e.g., ISchlegel et al.fll998|) and CO observations that 


of molecular gas (e.g., Dame et al.l I200H ). More recently, 


over the past few years, large-scale maps of diffuse interstel¬ 
lar bands tracing ( still unidentified ) large molecules have 
started to emerg e (lLan et al ] I20l4 IZasowski et al.l 120151 : 
iBaron et alj|2015f ). Surprisingly, the large-scale distribution 
of metals in the atomic gas phase - a fundamental compo¬ 
nent of the ISM - remains poorly measured. Characteriz¬ 
ing their spatial distribution is important: while metals are 


expected to represent about 2% of the mass of the inter¬ 
stellar medium (ISM), they often determine the chemistry, 
ionization state and temperature of the gas. Their spatial 
distribution reflects the return to the ISM of gas that has 
been processed in stars and stellar explosions. As of today, 
no large-scale map of metals in the gas phase exists. 


Gaseous metals can be probed using absorption line 
spectroscopy. Among the features accessible in the optical 
window, the Ca II and Na I lines are the most prominent. 
They correspond to the H, K and D lines discovered by 
Fraunhofer in the solar spectrum two hundred years ago 
(Fraunhofer 1814). Over the last two centuries, these Fraun¬ 
hofer lines have played an important role in the understand¬ 
ing of a varie ty of astronomica l phenomena both in the M ilky 
Way (e.g., iHartmarml Il904 
on extragalactic s cales (e.g., 


Routlv & Spitzcr 1952) and 


Bokscnbcrg & Sargent] 1 19781 : 


iBlades et al.lll98ll) . In this paper we use a large set of ex¬ 
tragalactic spectra observed by the Sloan Digital Sky Sur¬ 
vey (SPSS. [York et aklhood ) to measure the absorption in¬ 
duced by the Milky Way. Interestingly, while such informa¬ 
tion is encoded in virtually any extragalactic spectrum, no 
systematic extracti on and analysis h as been attempted on 
a large scale. Only IPoznanski et al.l d2012l ) measured Na I 
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absorption in SDSS galaxy spectra to characterize its cross¬ 
correlation with dust reddening. We create maps 0 of Ca II 
and Na I absorption over a quarter of the sky induced by 
gas in and around the Milky Way, i.e., in both the ISM and 
the circumgalactic medium (CGM), and study their global 
correlations with other baryonic tracers, such as neutral hy¬ 
drogen and dust grains. The Ca H & K doublet corresponds 
to the fine structure splitting of the singly ionized calcium 
excited states Ca II (Ca + ) at A = 3934.78 A (K) & 3969.59 A 
(H) in vacuum, and Na D, which is also a doublet, corre¬ 
sponds to the fine structure splitting of the neutral sodium 
excited states at A = 5891.58 A (D 2 ) & 5897.56 A (Di) in 
vacuum. Since their discovery they have played important 
roles in astrophysics because of their strength and their lo¬ 
cation in the visible part of the spectrum. Both calcium and 
sodium are among the most abundant heavy elements, with 
log(Ca/H)Q + 12 ~ 6.34 and log(Na/H)o + 12 ~ 6.24 (As- 
plund 2009), and are found in a variety of astrophysical en¬ 
vironments. The ionization potentials are 6.11 eV, 11.87 eV 
and 50.91 eV for Ca I, Ca II and Ca IH, and 5.14 eV and 
47.29 eV for Na I and Na II (lMortonll2003f ). With ionization 
potentials lower than 1 Rydberg, Ca II and Na I are therefore 
not the dominant stage of ionization in most astrophysical 
environments. 

The paper is organized as follows. In Section 0 we 
briefly describe the data set and method. In Section [3] we 
present the maps of Ca II and Na I absorption and inves¬ 
tigate the correlations of metal absorption with neutral hy¬ 
drogen and dust. We discuss the effect of metal absorption 
on precision photometry in Section [4] Section [5] summarizes 
our main results. 


2 THE DATA 


We derive our absorption measurements from quasar and 
galaxy spectra taken f rom the SDSS Data Release 7 (DR7, 
lAbazaiian et al.l 12009). For quasars, we u se the quasar 
catalog compiled by Schneider et ah d2010l) and t he im - 
proved redshift estimates given by Hewett fe Wildl J20ldl . 
For galaxies, we use th e MPA-JHU value-added catalog (e.g., 
iTremonti et al] l2004T l. Absorption line measurements re¬ 
quire accurate estimation of the intrinsic background source 
continuum F( A). Below we briefly describe the continuum 
estimates we use: 

• Quasar con tinuum estimation: We use the esti¬ 
mates provided bv lZhu fe Menardl d2013alh who used a basis 
set of quasar eigenspectra created with the dimensionality- 
reduction technique nonneqative matrix factor ization (NMF, 
iLee fe Seund 1 19991 : iBlanton fe Rowelsl l2007tl . Large-scale 
residuals not accounted for by the NMF basis set are re¬ 
moved with appropriate median filters. This set of flux 
residuals has been used to create a sample of about 
50,000 absorber systems dZhu fe Menardl l2013al b to mea¬ 
sure the_total_amount_of Ca II around low-redshift galax¬ 
ies dZhu fe Menardl 2013bj)_and Mg II around luminous red 
galaxies at z ~ 0.5 ( Zhu et ai1l2014 h 


• Galaxy continuum estimation: We use estimates 
of the ir intrinsic continuum fluxes provided by IZhu et all 
d20ld) . These authors modeled the observed galaxy spec- 
tra using single stellar population (SSP) models of 
iBruzual fe Charlotl d2003l l with the Padova 1994 library of 
stellar evolution tracks (e.g., lAlongi et aid fl99^1 and the 
IChabrierl d2003d) initial mass function (IMF). For this anal¬ 
ysis, we restrict the selection of galaxies to those at redshift 
> 0.2, most of which are so-called luminous red galaxies 
(LRGs) with well-understood spectral energy distribution. 
Because red galaxies have low-level fluxes blueward of the 
4000 A break, there is limited information of Ca II absorp¬ 
tion induced by the foreground sources, we will use galaxy 
spectra only for Na I absorption. 

In total, we make use of spectra of 84, 533 quasars and 
181, 286 galaxies. These sources cover about a quarter of 
the sky, mostly at high Galactic latitude in the northern 
hemisphere with relatively low dust extinction. We note 
that all the background sources used here are extragalactic 
objects and the spectral resolution is 69 km s -1 . All the 
absorption measurements presented here therefore include 
contributions from all gas, including both in the ISM and 
the CGM, of the Milky Way. 


To cross-correlate our maps of metal absorption with 
the distributions of hydrogen and dust we use the following 
datasets: 


• Neutral hydrogen: The Leiden/Argenti ne/Bonn 
(LAB) Galactic H I Survey iKalberla fe Kerpl (|2009l) merged 
the L eiden/Dwingeloo Survey fLDS: lHartmann fe Burton I 
Il997f) in the northern sky an d the Instituto Argentino de Ra- 
dioas tronomfa Survey fIAR: lArnal et aI1l200ol : iBaiaia et al.l 
l2005ll in the south and present an all-sky brightness temper¬ 
ature map of the 21 cm hyperfine emission by neutral hydro¬ 
gen. The LAB Survey is the most sensitive Milky Way H I 
survey to date, with the most extensive coverage both spa¬ 
tially and kinematically. In the optically thin regime, the in¬ 
tegrated emission yields total neutral hydrogen column den¬ 
sity Ah 1 in the Milky Way, which we use to compare with 
our statistical metal absorption measurements. For Ahi, we 
use the values integrated over the full velocity range, from 
—450 km s' 1 to 400 km s -1 . 

• Dust: [Schlegel et alj (1 19981 '! derived Galactic reddening 
Eb-v over the full s ky from a composi te 100/rm of the 
COBE/DIR BE (e.g.. lEtoggess et al.lll9~9l) and IRAS/ISSA 
maps (e.g., IWheelock et al. Ill994h . assuming uniform dust 
properties with values of Rv = Av / Eb-v « 3.1. We here 
use Eb-v as a proxy of line-of-sight dust column density 
and adopt their derivations. 


For both Ahi and Eb-v, we use the maps distributed on 
the NASA LAMBDA websitfl which were produced under 
the HEALPix scheme (iGorski et al. Il2005lh The HEALPix 
pixelization divides a spherical surface into pixels with the 
same area. The main parameter of the pixelization is the 
number of pixels along a fixed longitude -Aide. Both Ahi 0 
and Eb-v maps were stored with A s id e = 512, which divides 


2 http://lambda.gsfc.nasa.gov/product/foreground 

1 The maps of Ca+ and Na absorption are made publicly avail- 3 These data were used in lLand fe Slosarl d2007r) and given for 

able at: http://www.inasan.ru/~khramtsova/project_CaNa public use through LAMBDA. 
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Quasar spectra 


Galaxy spectra 




Figure 1. Examples of stacked continuum-normalized residual spectra of quasars (left panel) and galaxies (right panel) in four Galactic 
latitude bins. The red lines represent the best-fit double-Gaussian profiles. 


the whole sky into 3,145,728 pixels, with the resolution of 
approximately 0.01 sq. deg per pixel. 


3 ANALYSIS 

Given the typical signal-to-noise ratio of SDSS quasar and 
galaxy spectra, the Ca II and Na I absorption lines induced 
by the Milky Way are in general not detectable in individ¬ 
ual spectra. Their detection requires constructing averaged 
continuum-normalized residual spectra. To do so we use an 
inverse-variance weighted estimator, using the wavelength- 
dependent errors given by the SDSS pipeline. With the 
stacked residual spectra we then measure the rest equivalent 
width of the Ca II and Na I doublets with a double-Gaussian 
profile, allowing the width and line ratio to be free param¬ 
eters. In Figure Q] we present examples of stacked spectra 
as a function of Galactic latitude, showing clear detections 
of Ca II and Na I absorption. The best-fit double-Gaussian 
profiles are presented with the red lines. As shown in the 
right panel of the Figure, there is limited information about 
Ca II absorption using galaxy spectra. This is because we 
are using red galaxies which have little fluxes blueward of 
the 4000 A break. In the following, we therefore use only 
quasars as background sources for Ca II, while for Na I we 
combine both quasars and galaxies as a whole sample. 

3.1 Absorption Maps 

To create maps of Ca II and Na I absorption we pixelize 
the sky using the HEALPix scheme. We choose the num¬ 
ber of pixels along the Galactic longitude to be lV s ide = 16, 
which divides the whole sphere into 3072 pixels with equal 
area of about 3.7° x 3.7°. This choice of resolution is mo¬ 
tivated by the surface number density of objects so that in 
the majority of the pixels, we have enough objects for a 
signal-dominated determination of the absorption strength. 


In total, our data cover about 700 pixels (about 9,100 deg 2 ), 
and we have on average about 120 quasars and 270 galaxies 
in a pixel. For each of these resolution elements we perform a 
double-Gaussian fitting to the stacked residual spectra and 
measure the rest equivalent width of each line. At that res¬ 
olution, the typical S/N with which such lines are detected 
is about 3. 

Figure [2] presents maps of Ca II K and Na I D 2 ab¬ 
sorption in the Mollweide projection. They show the equiv¬ 
alent width of each line from 0 to 0.5 A using a 3-degree 
Gaussian smoothing. Note that to visually enhance the con- 
strast, equivalent width values greater than 0.5 A are satu¬ 
rated. This affects only a small fraction of the pixels. Both 
maps show large-scale gradients as a function of Galactic 
latitude. The distribution of Na I appears to be more con¬ 
centrated towards the disk of the Milky Way. To characterize 
this dependence, we measure the mean absorption strength 
as a function of latitude, using only data from the Northen 
Galactic hemisphere where the coverage is more complete. 
We present the mean equivalent width of Ca II and Na I 
measured in bins with A b = 10° in Figure [3] where the 
errorbars show the la dispersion of all the pixels. For com¬ 
parison, we also show the average H I column density in the 
bottom panel, where Ah i is the median value over all pixels 
of the Ah 1 HEALPix map in each latitude bin within the 
same SDSS footprint. Both the Ca II and Na I absorption, 
and Ah 1 decrease from low latitude towards high latitude. 
This is likely due to a combination of longer path through 
the disk at low latitude and higher local column density in 
the ISM. 

3.2 Correlations with Neutral Hydrogen and Dust 

More information on the distributions of Ca II and Na I 
can be extracted from these maps using cross-correlations 
with other tracers. This allows one to make use of signal 
lying below the noise level of individual HEALPix pixels. In 
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Figure 2. Maps of Ca II (top) and Na I (bottom) created from measurements of Fraunhofer H and D 2 absorption lines in the spectra of 
about 300,000 quasars and galaxies from the Sloan Digital Sky Survey. These maps use the Mollweide projection and show the equivalent 
width of each line from 0 to 0.5 A using a 3-degree Gaussian smoothing. Note that to visually enhance the constrast, equivalent width 
values greater than 0.5 A are saturated. This affects only a small fraction of the pixels. 


this section we present measurements of cross-correlations 
with the observed distributions of Galactic neutral hydrogen 
and dust. We note here that previous attempts to measure 
correlations between metal absorption, hydrogen and dust 
were typically based on datasets targetting “dense” clouds. 
By using spectra of extragalactic sources and integrating the 
absorption signal over a large number of lines-of-sight, our 
approach allows us to probe the total amount of absorption 
originating from both dense clouds as well as the diffuse 


interstellar and circumga lactic media. So f ar, on ly one such 
approach has been used: IPoznanski et al.l (120121 ) performed 
a similar statistical study of Ca II absorption. 

To investigate the dependence of Ca II and Na I ab¬ 
sorption on Ah i and Eg-v, we use the high-resolution 
HEALPix maps of Ahi and Eg-v distributed by the 
LAMBDA website. We divide all the HEALPix pixels into 
bins with Alog 10 A H i = 0.05 and Alog 10 A s _v = 0.05. 
We then find all objects (quasars and galaxies for Na I, 
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Parameters 

Ca K 

Ca H 

Na D 2 

Na Di 

Ao/10 21 (Ahi) 

13.18± 4.92 

32.53± 10.83 

2.08± 0.24 

2.36± 0.18 

«o (A H i) 

0.35± 0.03 

0.39 ± 0.03 

0.77 ± 0.04 

1.01 ± 0.04 

E 0 {Eg_v) 

4.21 ± 2.46 

10.90 ± 7.47 

0.39 ± 0.09 

0.26 ± 0.02 

ho ( Eg-v) 

0.29 ± 0.03 

0.33 ± 0.04 

0.63 ± 0.05 

1.06 ± 0.04 


Table 1. Best-fit parameters for the Wo~Ngi/Eg_y relations (Equation |T| and [2j 



Figure 3. Rest equivalent width of Ca II K (top panel) and 
Na I D 2 (middle panel) as a function of Galactic latitude b. The 
errorbars represent la errors from the double-Gaussian fitting. 
For comparison, in the bottom panel we show the average neutral 
hydrogen column density Ahi in the same footprint, using values 
from the LAB H I survey. The errorbars for Ah i are lc dispersion 
in the pixels used. 


quasars for Ca II) in the corresponding pixels, construct 
stacked continuum-normalized residual spectra and measure 
the equivalent width of Ca II and Na I. Note that we do not 
spatially downgrade the Ahi and Eg-v HEALPix maps 
(with iVside = 512) but instead we select objects in pix¬ 
els with the same Ahi or Eg-v values. In Figure [4) we 
present the correlations of the rest equivalent width of Ca II 
and Na I as a function of Ahi and Eg-v, where the error- 
bars are from the double-Gaussian fitting. Both Ca II and 
Na I absorption strengths correlate strongly with Ahi and 
Eg-v, increasing linearly towards higher density until satu¬ 
ration effect becomes important at Ahi~ 5 x 10 20 cm~ 2 and 


Eg-v ~ 0.08 mag. However the observed correlations indi¬ 
cate that the environme ntal dependence of Na I and Ca II 
abundance are different. iPoznanski et aD d2012l ) performed 
a similar correlation analysis between Na I and Eg-v and 
obtained similar result 

We characterize the trends by the following relations: 

w(Nm> = (itr)” A ■ m 

and 

W(E B -v)= A. (2) 

We estimate the best fit parameters considering only the 
range where Ahi< 5 x 10 20 cm -2 and Eg-v< 0.08. Above 
these values the observed trends indicate that the absorp¬ 
tion lines become saturated. The best-fit parameters are pre¬ 
sented in Tableland the corresponding trends are overplot¬ 
ted with dash lines in Figure [4] Equations [T] and [2] provide 
us with a way to estimate the average absorption strength 
of Ca II and Na I at a given Ah i and Eb-v- It is interesting 
to note that the slope of Wo -Ahi (and Eg-v) relation 
is close to unity, while the slope of Wo^H-Ahi is shallower 
with a best-fit value of about 0.4. The weaker dependence 
of the Ca II absorption on Ahi could be due to either the 
lower dust depletion level of calcium or higher ionization 
parameter or both in lower-density environments. We fur¬ 
ther investigate the abundances of Ca II and Na I and their 
dependence on Ahi in the next section. 

3.2.1 Abundance of Ca II and Na I 

We now investigate the dependence of Ca II and Na I abun¬ 
dances on hydrogen column density. We first estimate the 
Ca II and Na I column density using the weaker line of each 
doublet under the assumption of negligible saturation. Using 
the linear curve-of-growth relation at low column density, we 
convert the rest equivalent width of the weaker lines (W^ 1 
and Wq 1 ) to column density and calculate the relative abun¬ 
dance Acaii/AHi and ANai/AHi. When the absorption is in 
the linear regime of the curve of growth, the rest equivalent 
width is related to the column density with 

A = 1.13 x 10 20 cm" 2 (3) 

J* 

4 For reddening v alues Eb—v greater th an 0.2 mag, our results 
differ from those of lPoznanski et all (120121 ). We realized that these 
authors included measurements at high-extinction regions in the 
Galactic disk due to a few thousand stars erroneously included in 
the MPA-JHU “galaxy” catalog. Treating these sources as galax¬ 
ies contaminates the measurements at high E(B — V) values. 
These stars have been carefully removed from our analysis. 
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Figure 4. The dependence of Ca II and Na I absorption on neutral 
(Eb-v, right panel). The lines represent the best-fit linear relation* 



E(B-V) (mag) 


hydrogen column density (A^hi, left panel) and Galactic reddening 
with data at TVhi < 5 x 10 20 cm -2 and Eb—v < 0.08 mag. 




Figure 5. The dependence of Na I and Ca II abundance on neutral hydrogen TVhi- The open circles represent the measurements in our 
statistical analysis, integrated over the full velocity range. The dots and crosses show th e abundances of individual systems, HVC/IVC 
and ISM clouds, from high-resolution spectroscopy, compiled bv lWakker Sz Mat hid (|200(J) . The solid lines are their best-fit li near relations 
to the se individual systems. The dashed line in the left panel shows the best-fit relation to extragalactic Ca II absorbers by I Richter et all 
<201 if) . 


where both Wo and A are in unit of A and the os- 
cillator strength f for the Ca II H line is 0.324 
( Safronova ii ^ i SafronovaJbOllJ) and for Na I Di is 0.320 
( Kelleher fe Podobedoval 120081 ). This relation is valid when 
the optical depth is smaller than unity, which corresponds to 
iVcaii.Nai 10 12 ' 5 cm -2 (assuming the Doppler broaden¬ 
ing factor b ~ 10 km s _1 ). The measurements are presented 
in Figure[5]with open circles. The Na I abundance is roughly 
constant within the Nh i range probed, while the Ca II abun¬ 
dance decreases rapidly with increasing Nni- 

Our measurements are integrated values over the full 
velocity range with contributions from both the ISM in the 


disk and the CGM in the halo. To better understand the 
physics of the correlations, we compare our measurements 
with a compilation of abundances of i ndividual systems 
based on high-resolution spectroscopy dWakker fe Mathisl 
l200di . In Figure El the cross symbols represent the mea¬ 
surements of abundance in the ISM clouds, the dotted 
symbols show the measurements of abundance in high- 
velocity clouds (HVCs) and intermediate-velocity clouds 
(I VCs), and the solid lines are the best-fit linear relations 
by IWakker fe Mathisl 120001 1 for these individual systems. 
In the left pa nel, we also overplot the best-fit relation by 
iRichter et all 1201 ill who studied the abundance of Ca II 
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N(HI) (cm" 2 ) E(B - V) (mag) < N(Na I) > (cm- 2 ) 


Figure 6. The dependence of the N^a,!/NCa.ll ratio on JVhi, Eg—v and JVNai- 


of about ten extragalactic absorber systems. Our measure¬ 
ments of the average Na I abundance at a given total JVhi 
coincide with the abundance of individual systems, while the 
average Ca II abundance is closer to the relation of extra- 
galactic systems and overall about half a dex higher than 
those of ISM clouds, though with a similar dependence on 
JVhi. 

We would like to stress that the variables in our anal¬ 
ysis, the neutral hydrogen column density (JVhi), the metal 
absorption strength and column density, are all integrated 
values over the full velocity range. The neutral hydrogen col¬ 
umn density is dominated by the contribution from the ISM, 
because along a given line of sight the contribution from the 
ISM is orders-of-magnitude higher than from the CGM, even 
if the sightline covers an HVC or IVC. The Na I abundance 
of individual systems, though with a large scatter, on aver¬ 
age is roughly constant at different densities. The integrated 
Na I absorption in our analysis therefore must also be dom¬ 
inated by the contribution from the ISM and this explains 
why the average Na I abundance is consistent with that of 
the individual systems. The Ca II abundance of individual 
systems, on the other hand, is ord ers-of-magnitud e higher 
at low density than at high density. ]Zhu fe Menardl (|2013bl ) 
show that the total amount of Ca II in the CGM, averaged 
over all galaxies, is about an order-of-magnitu de larger than 
that in the ISM of the Milky Way (see also iRichter et al.l 
l201lll . In our analysis, making use of extragalactic sources, 
Ca II column densities are therefore expected to carry a 
significant CGM contribution, while hydrogen column den¬ 
sities are mostly originating from the Milky Way disk. We 
thus expect our measurements to display average values of 
Ca II abundance that are higher than those obtained for 
individual clouds, as shown in figure [5] 

In Figure [SJ we further investigate the dependence of 
the relative abundance of Na I and Ca II (JVNai/JVcaii) on 
Nh 1 , E B -v and JV N ai- The JVNai/JVcaii ratio is higher at 
higher hydrogen/dust/Na I density, increasing from about 
0.4 at JVhi ~ 10 20 cm- 2 ( Eg-v ~ 0.01 mag) to 2 at IVhi ~ 
10 21 cm -2 (Eg _y ~ 0.1 mag). We fit these relations with 
similar power-law functions as in Figure [U 


-NNal/Acall = 

/JV H i\“ 1 

V Ni ) ’ 

(4) 

and 



ANal/Acall = 

(^r • 

(5) 

and 




Parameters 

JVn a I/JVca II 

JVi/10 21 (JVhi) 

2.13± 0.40 

«i (JV H i) 

0.53± 0.04 

El ( Eg-v) 

0.22 ± 0.06 

hi ( Eg_v ) 

0.60 ± 0.07 

JV 2 /IO 12 (JV Na i) 

7.07± 0.82 

02 (JVNai) 

0.58± 0.03 


Table 2. Best-fit parameters for the A^ai/Ncaii ra tio- 
-^Hi/Fs—v/^Nai relations (Equation [4] [5] and |6|). 


AWiVcall = (^) 


( 6 ) 


The best-fit relations are given in Tableland overplotted in 
Figure [S] The slopes are both positive, indicating that the 
jV Na I /iVcja 11 ratio is higher in denser environments. 

This environmental dependence of the JVNai /Nca ii ra- 
tio has been noticed for decades (e.g., Phillips et ai .1119841 ; 
lYallerga et al.1 [l993l : iBen Bekhti et alj|2012h . However, the 
physical reason is yet to be determined. It is generally at¬ 
tributed to the different dust depletion levels of sodium 
and calcium at different densities, but other factors such as 
ionization effects can yield the same dependence. We here 

briefly discuss two theories. _ 

_ A possibl e scenario w as proposed by IPhillips et al.1 

(1 19841 . see also lBarlowlll978h : During the formation of dust 
grains, sodium is not highly depleted to begin with because 
of relatively high condensation temperature and low absorb¬ 
ing energy, so density does not play a role in its low deple¬ 
tion level. On the other hand, calcium can be easily accreted 
onto surface of dust grains due to relatively low condensa¬ 
tion temperature and high absorbing energy. At higher den¬ 
sity, its accretion rate is higher yet the destruction rate is 
lower because of greater UV shielding and enhanced man¬ 
tle growth, while at lower density, calcium trapped in grain 
surface is easier to be ejected by shock sputtering. All these 
dynamic processes result in higher calcium depletion level 
at higher density. 

When an analysis is sensitive to both ISM and CGM 
absorption, it is important to take into account of the envi¬ 
ronmental dependence of metallicity and dust abundance. 
Gas clouds in the CGM, such as HVCs and the Magel- 
lanic Stream ; have lower metallicity (about 0.1 solar , see 
IWakker et alJll^ fiRichto~_eta i]l2nnil : lFox et a.l .11201 .'ft ) and 

dust-to-gas ratio ( Lu et al.ll 1998h than those in the ISM. AI- 






















































8 Murga et al. 


though the total abundance of Ca and Na is lower in the 
CGM due to lower metallicity, the overall dust depletion 
level is also lower because of the low dust content. The com¬ 
bined effects of these dependences therefore can result in 
the observed dependence of the ANai/Ncaii ratio on the 
environment. 

Another theory is based on that the ionization pa¬ 
rameter is higher at lower density, and neither Ca II nor 
Na I is the dominant ionization state. Models with density- 
dependent ionization parameters but constant dust deple¬ 
tion level can reproduce the dependence of the ANai/Vcaii 
ratio on environment 0. The exact mechanism responsible 
for this density dependence of Ca II and Na I abundances is 
still an open question. 


3.3 The Routly-Spitzer effect measured on large 
scales 

iRoutlv fe Spitzerl (jl952) found that the average AN ai /Nca. ii 
ratio decreases when the cloud velocity |ulsr| (measured in 
the local standard of rest) increases. This so-called Routly- 
Spitzer effect is usually derived from measurements of indi¬ 
vidual s ystems in the local IS M within te ns of kilometer per 
second (IVallerga et al.l 1 19931 1. Recently, iBen Bekhti et ahl 
(l2012ll reported a similar effect for individual absorption sys¬ 
tems, based on deviation velocity measurement^] in the halo 
of the Milky Way, extending the detection of the Routly- 
Spitzer effect to larger scales. Here we investigate if such an 
effect also exists in our statistical data. 

To study the dependence of the ANai/AIcaii ratio on 
i>lsr, we use the HEALPix maps as presented in Figure[2] In 
each HEALPix pixel, we determine the centroid wavelength 
from the double-Gaussian fitting and calculate the veloc¬ 
ity offset from the rest-frame wavelength of the doublets. 
When doing so we find that Na I and Ca II are on average 
blueshifted by 8 and 23 km s _1 in the heliocentric reference 
frame. We then convert these observed velocities to the loca l 
standard of rest (LSR, ICourteau fc van den Bergh Ill999l h 
We show the corresponding distributions of absorption ve¬ 
locities for Ca II and Na I in the northern hemisphere in 
Figure [3 As can be seen, Ca II tends to be found at lower 
velocities compared to Na I. We find the mean of the Na I 
velocity distribution to be about —3 km s^ 1 , while the value 
for Ca II is found to be about —17km s -1 . Recent H I sur¬ 
veys have shown that, in the northern sky, the high- and 
inter mediate- velocity clouds have on average n egative veloc¬ 
ities dRichter et al.l [20031 : IPutman et aTl 1 201 2l 1 . The higher 
detection rate of Ca II at negative velocities is therefore 
consistent with the interpretation that Ca II detected in the 
extragalactic sources traces a significant amount of gas from 
the CGM. 

We then group the pixels with similar Vlsr values and 
measure the mean W / o' aI and Wo CaI1 in each Ulsr bin. We 
note that using the velocity derived from Ca II produces 


5 Bart Wakker, private communication 

6 Deviation velocity is the difference of the radial veloci ty of the 
clou d and the terminal velocity of the Milky Way disc dWakkeri 
Il99lh . However, for high galactic latitude systems, the effect 
should be similar regardless of the rest frame. 



Figure 7. The velocity distributions of Ca II and Na I absorp¬ 
tion in the local standard of rest using measurements from the 
northern Galactic hemisphere. The distribution for Ca II is offset 
from the origin and wider than that of Na I. This is consistent 
with the interpretation that Na I originates mostly from the disk 
while Ca II also traces a significant amount of material from the 
circumgalactic medium. 


similar trends with a larger scatter. The left panel of Fig¬ 
ure [8] shows our measurements of the ANai/ANaii ratio as 
a function of Ulsr- Interestingly, we find that the average 
ANai/ANaii ratio shows a similar dependence on vlsr as the 
local Routly-Spitzer effect for individual systems, therefore 
our data show that the velocity depe ndence of ANai/Acai i 
ratio also exists on the Galactic scale. iBarlow fe Silkl (Il977l ~) 
suggested that it is related to dust destruction. One of the 
major mechanisms of dust destruction is sputtering by in¬ 
terstellar shocks. The destruction rate depends on the shock 
velocity which results in higher dust destruction rate and 
higher Ca II abundance relative to Na I at higher velocity. 
Subsequent theoretical analyses and simulations have lent 
support to su ch proposition (e.g ., iDraine fe Salpeter~lll979l : 
iMcKeel 1 19891 ; I Jones et al.l 1 1994 ) . However, it is also possi¬ 
ble that this effect can be just a consequence of the density 
dependence of the ANai/IVcaii ratio shown in the previous 
section, regardless of the actual cause. In the middle and 
right panels of Figure [5] we show that the Na I density 
has a similar dependence on vlsr as the ANai/ANaii ratio, 
while Ca II density has a much weaker dependence. Since the 
ANai/ANaii ratio is higher in higher-density environment, as 
probed by Na I, H I and dust (Figure [5]), this can therefore 
explain part of the velocity dependence. The mechanisms 
responsible for the Routly-Spitzer effect on large scales can 
therefore be more complicated than interstellar shock only. 


4 IMPLICATIONS FOR PHOTOMETRIC 
OBSERVATIONS AND REDSHIFT 
ESTIMATION 

The Ca II AA3934, 3969 and Na I AA5890, 5896 absorption 
lines due to the ISM and CGM of the Milky Way are ex- 
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Figure 8. The global Routly-Spitzer effect. The left panel shows the dependence of the i /Nc\- A \[ ratio on dlsr of Na I. The right 
panels show IVNal arid Nq & h as a function of ulsr- Note the y-axis ranges in the right two panels are consistent, i.e., the maximum is 
a factor of four of the minimum in both cases. 


pected to impact the broad band photometry of extragalac- 
tic objects. Two types of effects can occur: 

• absorption on top of continuum flux: Metal line absorp¬ 
tion is expected to suppress a fraction of the expected flux 
from extragalactic sources. This effect has been statistically 
detected in the broad band photometry of distant quasars 
iMenard fe Fukugita I d2012l 'l. We note that the distribution 
of metals in the gas phase does not exactly track that of the 
dust. Therefore infrared emission (as used in SFD) cannot 
be used to precisely map this kind of absorption. 

• Absorption at the location of emission lines: if a red- 
shifted emission line happens to coincide with one of the 
Galactic absorption features, the apparent broad-band color 
of the corresponding galaxy can be appreciably modified. 
This can happen when, for example, [O III] A5007 is seen on 
top of Na D at z ~ 0.176 or [O II] A3727 is seen on top of Na 
D at 2 ~ 0.58. This introduces a level of degeneracy between 
redshift and color for star forming galaxies and is expected 
to impact the accuracy of photometric redshift estimation 
(see Rahman et al. (in prep) for a more extended discussion 
of this effect). 


5 SUMMARY 


• We have mapped out the Na I D and Ca II H & K 
absorption over about a quarter of the sky, obtaining signal- 
dominated measurements at a resolution of about one de¬ 
gree. The corresponding maps are publicly available. 

• Cross-correlating these absorption maps with H I col¬ 
umn density and dust reddening maps, we find that the de¬ 
pendence of Ca II absorption on Ahi and Eb-v is much 
weaker than Na land provide fitting formulae for the average 
relations between these quantities. We find the abundance 
of Na I with respect to neutral hydrogen to be roughly con¬ 
stant in different environments, while the Ca II abundance 
decreases with hydrogen column density. We discuss the pos¬ 
sible mechanisms that can cause such different behaviors, 
including dust depletion and ionization effects. 

• We show that, on average, the IVNai/Acaii ratio de¬ 
creases with velocity with respect to the local standard of 
rest up to velocities reaching 100 km s -1 . Interestingly, this 
is consistent with the local Routly-Spitzer effect originally 
measured for individual systems within a much smaller ve¬ 
locity range. 

• We commented on the fact that Galactic absorption 
lines can change the apparent colors of emission line galaxies 
at specific redshifts, which can possibly impact photometric 
redshift estimation. 


The light emitted by extragalactic sources has to penetrate 
the ubiquitous gas and dust in the interstellar and circum- 
galactic media of the Milky Way. Every extragalactic spec¬ 
trum therefore carries information about these two environ¬ 
ments. 

We have used the large collection of quasar and galaxy 
spectra obtained by the Sloan Digital Sky Survey to detect 
absorption lines induced by the presence of Ca II and Na I. 
The absorption signal is typically too weak to be detected in 
individual SDSS spectra. Measuring it requires a statistical 
appro ach. Not until recently had this avenue been explored 


& 


Poznanskietal . 20121 : lLan et al.l[2014l : IZasowski et al.l 


Baron et al.l 120151 1. We have created maps revealing 


the spatial distribution of Ca II and Na I absorption on the 
sky. We have then used these new datasets to study corre¬ 
lations between these two components and other baryonic 
tracers. Our analysis has allowed us to obtain the following 
results: 


The detection of absorption features relies on the ability 
to accurately estimate the intrinsic spectral energy distribu¬ 
tion of the background source. To do so w e have used the 
continu um estimates of quasars presented in lZhu fe Menardl 
(l2013al l. b ased on a vector- decomposition method, and of 
galaxies bv IZhu et al.l J20ldi 'l. based on stellar spectral tem¬ 
plates. Our work demonstrates that large spectroscopic sur¬ 
veys, though designed for studies of the light sources them¬ 
selves or cosmology, provide valuable datasets for studies of 
the ISM and CGM of our own Milky Way. We have only just 
started to explore this avenue and there is more information 
to be extracted. For example, the absorption induced in the 
spectra of extragalactic sources includes contributions from 
both the ISM and CGM of the Milky Way, while absorption 
in the stellar spectra only has contribution from the ISM 
up to the star’s distance, differential analysis therefore can 
be performed to produce a 3D map of the gas distribution 
in and around the Milky Way. The continuum-estimation 
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methods and statistical absorption-line spectroscopic ap¬ 
proaches we developed are robust, generic, and readily appli¬ 
cable to any large datasets from ongoing and future surveys, 
such as SDSS-IV, DESi and PFS. 
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